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g : ABSTRACT 

(N 

rH ■ We present the first report on chaos synchronization between two unidirectionally linearly and non- 

^ . linearly coupled systems with multiple time-delays and using the Razumikhin-Lyapunov approach 
find the existence and stability conditions for different synchronization regimes. The approach is 
tested on the famous nonlinear models-Ikeda, Mackey-Glass and Lang-Kobayashi systems. 
PACS number(s):05.45.Xt, 05.45. Vx, 42.55.Px, 42.65. Sf 
Q ■ Introduction- There is continued growth in the field of chaos control [1] and ever-increasing 

rj . appreciation of its appications among researchers. Application of chaos control can be found in 
secure communication, optimization of nonlinear system performance, modeling brain activity and 
pattern recognition phenomena [2], species population control [3], etc. 

Recently delay differential equations (DDEs) have attracted a lot of attention in the field of 
nonlinear dynamics. DDEs are used to model dynamical systems in many scientific and enginering 
^ ■ areas,e.g. optics, biology, climatology, economy, cryptosystems based on synchronized hyperchaos, 
O ■ networks [4-5], etc. In comparison with a single time-delay DDEs with multiple time-delays are 
^ \ more relaistic models in the interacting complex systems. These delays result naturally from 
\Q \ the finite propagation velocity of information, from the latency of feedback loops, from the finite 
O ■ switching times between different states of the system. Additional time-delays could be useful 
(^ ■ e.g. to stabilize nonlinear system's output [6]. The application possibilities based on chaos re- 
quire proper control of complexity. To the best our knowledge chaos synchronization between 
the systems with multiple time-delays has not been investigated yet. Having in mind enormous 
n \ application implications of chaos synchronization e.g. in secure communication, performance 
,^ ' optimization in nonlinear systems, stabilization problems, etc. investigation of synchronization 
^ ■ regimes in the multiple time-delay systems is of paramount importance. 

In this Letter we present the first report on chaos synchronization between two unidirec- 
tionally linearly and nonlinearly coupled chaotic systems with mutiple time-delays and find the 
existence and stability conditions for different synchronization regimes. We test the approach on 
the paradigm Ikeda, Mackey-Glass and Lang-Kobayashi models [4]. We hope this research will 
pave the way for the intensive experimental investigations of chaos synchronization in the multple 
time-delays systems. 

General approach-Consider synchronization between the double-feedback systems of general 
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form , 

dx 

= -ax + mif{xrj + m2f{xr^), (1) 

dv 

— = -ay + m^fiyrj + mJiy^^) + Kf{xr^), (2) 

where / is differentiable generic nonlinear function. Throughout this paper Xr = x{t — r). 
One finds that under the condition 

rrii — K = 1713,1712 = m4^ (3) 

Eqs. (1) and (2) admit the synchronization manifold 

y X-j-^—r-^ . y^J 

This follows from the dynamics of the error A = x^g-ri — y 

dA 

— = -aA + Tn^ArJ'iXr,,) + m2ArJ'{Xr.,+rz-T^)- (5) 

Here /' stands for the derivative of / with respect to time and the derivative should be bounded. 
The sufficient stability condition of the trivial solition A = of (5) can be found using Razumikhin 
theorems (see for details, [5], pp. 151-161). In [5] considering the Lyapunov function V = ^A^ and 
using Razumikhin theorems (hereafter we use the term Razumikhin-Lyapunov approach), under 
restrictive inital data condition |A(t)| > |A(t-|-^)|, for all 6 G [—Tj, 0] it is shown that zero solition 
of ^ = — a(t)A(t) — Y^'jZibj{t)A{t — Tj{t)) is uniformly asymptotically stable for all bounded 
continous functions a,bj,Tj if a{t) > 6 > 0,J2'jZibj{t) < p6,0 < p < 1,0 < Tj{t) < r for all 
t G (—00,00). Thus, applying the Razumikhin-Lyapunov approach we obtain that the sufficient 
stability condition for the synchronization manifold (4) can be written as: 

a > |m3(sup/'(x^3))| + |m2(sup/'(a:^2+^3_^J)|. (6) 

Here sup f'{x) stands for the supremum of /' with respect to the appropriate norm. 
We notice that for T3 > ti,T3 = ti, and T3 < ti (4) is the retarded, complete, and anticipating syn- 
chronization manifold [7,8], respectively. Analogously one finds both the existence (m2 — K = 777,4, 
"7i = ms) and sufficient stability (a > |7?73(sup /'(x^-g))! + |7?72(sup /'(a:T-2+r3-Ti))|) conditions for 
synchronization manifold y = Xr^^^r2- 

In the case of linear coupling of form K{x — y) one obtains that under the condition mi = 
^3, 'm-2 = 'nii synchronization manifold y = x exists and it is stable ii a + K > | mi (sup /'(xtJ) | + 
|m2(sup /'(xt-j))]. For the parameter mismatches, e.g. ri 7^ ri/, (ri/ is the delay time for the first 



feedback loop in y dynamics) it is clear that complete synchronization is not the synchroniza- 
tion manifold. Then for such a case we use the auxiliary system method to detect generalized 
synchronization [9]:that is given another identical driven auxiliary system z{t), generalized syn- 
chronization between x{t) and y{t) is established with the achievement of complete synchronization 
between y{t) and z{t). Thus, the auxiliary method allows to find the local stability condition of 
the generalized synchronization [9]. Applying this method we find local stability condition of the 
generalized synchronization between y and x\ a + K > |m3(sup /'(?/t-i/))| + |"^4(sup /'(|/t-2))|- We 
also notice for Kxr^ type of coupling with mismatches between the relaxation coefficents the only 
possible synchronization manifold is y = Xr^ (independent of the relation between the coupling and 
feedback delay times) with existence a2 — ai = K,mi = nis, and m2 = m^ and sufficent stability 
conditions 0:2 > | mi (sup /'(a^n+rs))! + |m,2(sup /'(a:r2+r3))|- The presence of such a mismatch could 
be useful in the interpretation of the future experiments with coupling-delay lag synchronization 
and could serve as a "switch off' mechanizm for certain types of synchronization manifolds. 
Generalization of the approach to n-tuple feedback systems, i.e. systems with multiple delays 
of type (1) and (2) is straightforward. We underline that a stability condition derived from the 
Razumikhin-Lyapunov approach is a sufficient condition: it assures a high quality synchronization 
for a coupling strength estimated from the stability condition, but does not forbid the possibility 
of synchronization with smaller coupling strengths. The threshold coupling strength can be esti- 
mated by the dependence of the maximal Lyapunov exponent A of the error dynamics on K:i.e. 
from X{K) = [10]. 

Example l:The Ikeda model- First we test the approach on the nonlinearly coupled Ikeda 
model: ^ = —ax + m,i sin x^j + m,2 sin x^-j ; ^ = —ay + m^, sin y^-^ + 7714 sin y^.^ + K sin Xr^ , with pos- 
itive ai^2 and — mi^2,3,4-This investigation is of considerable practical importance, as the equations 
of the class B lasers with feedback (typical representatives of class B are solid-state, semiconductor, 
and low pressure CO2 lasers [11]) can be reduced to an equation of the Ikeda type [12]. The Ikeda 
model was introduced to describe the dynamics of an optical bistable resonator, plays an important 
role in electronics and physiological studies and is well-known for delay-induced chaotic behav- 
ior [7,13]. We find that the Ikeda systems can be synchronized on the synchronization manifold 
y = Xt-3-ti under the condition ttt-i — K = rris, m2 = m^ and is stable if a > [msl -|- |r?7,2|. Numerical 
simulations fully support the analytical results. The Ikeda model was simulated using the DDE23 
program [14] in MATLAB 6. Figure 1 shows the time series of the driver x{t) (solid line) and 
the driven system |/(t) (dotted line) for a = 5, ri = 3,r2 = 2,r3 = 1, m2 = m4 = —1,7711 = —18, 
m.3 = — 1 and K = —17. After transients the driven system shifted T3 — T1 = —2 time units to the 
left and y = x{t + 2) (anticipating synchronization). 

Example 2: The Mackey- Glass model.-ConsideT complete synchronization between the lin- 
early coupled double- feedback Mackey-Glass systems: ^ = — aix -|- ki^^ — h k2^T^; ^ = 



^2 



— «2?/ + ^3Tx^s — \- ki-r^ — h K[x — |/).The dynamical variable in the Mackey-Glass model is the 
concentration of the mature cells in blood at time t and the delay time is the time between the 



initiation of cellular production in the bone marrow and the release of mature cells into the blood 
[1,10]. (At present there is also an electronic analog of the Mackey-Glass system [10].) We find 
that the Mackey-Glass systems can be synchronized on the synchronization manifold y = x under 
the existence ki = k^, k2 = k^ and stability a + K > {ki + ^2) ^^^ conditions. For analytical 
estimations we took into account that the absolute maximum of the function \f'{xr)\ is obtained 
at Xr = (^ry)^ and is equal to ^ ~^^ [10]. Figure 2 shows numerical simulation of the linearly cou- 
pled Mackey-Glass models:time series of the driver x{t) (solid line) and driven system y[t) (dotted 
line). The parameters are ti = 14, T2 = 20, a = l,b = 10,ki = k^ = 2, k2 = k4 = 0.2, K = 5. After 
transients the driven systems trajectory completely coincides with that of the driver system. In 
Figure 3 generalized synchronization between the linearly coupled Mackey-Glass models is shown 
for Ti = 14, Tif = 16;the other parameters are as in Fig.2. 

We emphasize that as the coupling strength estimated from the stability condition gives a high- 
quality synchronization, the synchronization manifold is robust against perturbations of the cou- 
pling strength. As mentioned above the onset of synchronization occurs at the coupling stength 
when the maximal Lyapunov exponent of the error dynamics vanishes as function of K [10]. Our 
estimations show that for the parameters values as in Fig.2 the threshold value of K is:i^ ~ — 1.11, 
which is far less than K = 3.44 found from the sufficient stability condition. 

Example 3:The Lang-Kobayashi model- As the last example we study synchronization between 
the unidirectionally coupled double-feedback Lang-Kobayashi systems [1,6-8] (laser diodes with 
a double external cavities, see [6] and references there-in): —^ = - — y^(Gi,2(A'^i,2 — ^01,02) — 
7i,2)^i,2 + A;i,2^i,2(^-T-i)exp(-/,u;ri)-h/i;3,4^i,2(t-r2)exp(-/,u;r2)-F-/i;5^i(t-r3)exp(-iu;r3);^^ = 
•h,2 — 7ei,e2A^i,2 — Gi,2{Ni^2 — A^oi,02) |-E'i,2|^) where £'1^2 are slowly varying complex fields for the 
master and slave lasers, respectively A'^1^2 are the carrier densities;A^oi,02 are the carrier densities at 
transparency;7i,2 are the cavity losses; ai,2 are the linewidth enhancement factors; Gi^2 are the op- 
tical gains; ki^2 and A;3^4 are the feedback levels for the master and slave lasers, respectively. /^s is the 
coupling rate;a; is the optical feedback frequency without feedback; ri,2 are the round-trip times in 
the external cavities for the coupled lasers; r^, is the time fiight between the master laser and slave 
laser (coupling delay time);Ji 2 are the injection currents; 7^i^e2 are the carrier lifetimes. The term 
/c5 exists only for the slave laser. In order to find possible synchronization regimes we compare 
e.g. equations for the dynamics of E2 and A''2 with dynamics of E'l^T-s-n and A^^i^T-a-n and find that 
for the case of identical lasers (all parameters are the same except for the feedback and coupling 
rates) Lang-Kobayashi systems can be synchronized on the synchronization manifold I2 = Ii^rs-n 
if 

ki = k3 + k5,k2 = k4^ (7) 

as the intensities I2 and Ii^r-^^n {I = \E\'^) can be made identical under these conditions. We have 
found the same existence conditions for the unidirectionally coupled laser diodes with incoherent 
optical feedbacks [15]. Unfortunately the stability of the synchronization manifolds for the Lang- 
Kobayashi model practically can not be studied analytically, and therefore one have to rely on the 
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numerical methods. Figure 4 shows synchronization manifold I2 vs./i^T-3-ri for parameters values 
A^o = 1-7 X 10^G' = 2.14 X 10^ n = 10~^s,r2 = 1.5 x 10~^s,ts = 2 x 10^^ a = 5, 7e = 0.9 x lO'^s, 
^ = 635nm,J = 0.02-fNo, ki = lOns-^h = lns-\k5 = 9ns^\k2 = ki = lOOns"!. 

Conclusions- By using the Razumikhin-Lyapunov approach we have presented the first report 
on different synchronization regimes between two unidirectionally coupled (linearly and nonlin- 
early) chaotic systems with multiple delays. We have successfully applied the approach to the 
paradigm models in nonlinear physics-the Ikeda, Mackey-Glass and Lang-Kobayashi models. We 
have found analytically the existence and whenever possible stability conditions for the anticpat- 
ing,lag, complete and generalized synchronization regimes. We hope that this research opens up 
possiblities for highly anticipated intensive experimental investigations of chaos synchronization 
in systems with multiple time-delays. 
The author acknowledges support from the Abdus Salam ICTP (Trieste, Italy) Associate scheme. 

Figure captions 

FIG. 1. Numerical simulation of the Ikeda model: the time series of the driver x{t) (solid line) and 

the driven system ?/(t)(dotted hue);?/ = x{t + t^ — ti) = x{t + 2). Dimensionless units. 

FIG. 2. Numerical simulation of the Mackey-Glass model: complete synchronization between y and 

X. Dimensionless units. 

FIG. 3. Numerical simulation of the Mackey-Glass model: generalized synchronization between y 

and X. Dimensionless units. 

FIG. 4. Numerical simulation of the Lang-Kobayashi modehthe dependence of I2 on Ii^-^^^^^. 
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